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Abstract. This article presents a statistical-mechanical 
treatment of a relationship (the force equation) between 
the gravitational potential for two particles and the corre¬ 
lation functions in a relaxed distribution of particles with 
different masses. This relationship is used in the case of 
galaxies interacting through a Newtonian potential in an 
Universe in expansion, i.e. the large-scale distribution of 
galaxies. 

By applying this equation and from the observed two- 
point correlation function for galaxies as a —1.8 exponent 
power law, I derive the approximate dependence of a mass- 
mass correlation function as a —2.8 exponent power law, 
i.e. I infer that mass is more correlated than galaxies at 
short distances, when the distribution is considered as re¬ 
laxed. 


1. Introduction 

In a physical system consisting of a collection of particles, 
information can be obtained from the way in which the 
particles are distributed in space. For a Poissonian distri¬ 
bution in equilibrium there is no correlation among the 
the particles, the position of each one being independent 
on the position of the others. This distribution sets the 
particles randomly in space with equal probability for all 
the positions. In contrast, when the distribution deviates 
from Poissonian, the correlations among particles are not 
null; instead, the particles will be mutually interacting. 

The relationship between the correlations of the distri¬ 
bution and the interaction force is called the force equa¬ 
tion and can be obtained under a certain hypothesis (also 
called the cosmic virial theorem in the literature relating 
to large-scale structure; Peebles 1976). The development 
of statistical-mechanical theory of liquids has allowed this 
relationship to be achieved in the specific case where the 
particles are all equal, point-like and in Boltzmann equilib¬ 
rium, which is quite common for classical particles in ther¬ 
modynamic equilibrium (see for example Goodstein 1975, 


March & Tosi 1976). In general, statistical-mechanical 
tools are applied to obtain the distribution of particles 
(molecules on the atomic scale) from an assumed form 
of the interaction (the Lennard-Jones interaction among 
molecules). 

Since these statistical applications can be developed in 
environments other than on the atomic scale, my purpose 
is to derive a few relationships that will be useful in cos¬ 
mology. Groups of galaxies in the large-scale structure of 
the Universe are observed in sky surveys, i.e. we have a 
distribution of particles in space although there are some 
differences with respect to the liquid. First of all, we know 
that these particles are not all identical; their mass and 
other characteristics that are not important dynamically 
are distinctive among themselves. In view of this, we could 
face problems when applying statistical mechanics on as¬ 
tronomical scales, but the issue of different masses is still 
an important difficulty to resolve hrst. 

Our purpose is to obtain a force equation for a dis¬ 
tribution of particles with different masses that will allow 
different magnitudes of interest to be related among them¬ 
selves. Our result will apply directly to Newtonian gravity 
for a Universe in expansion. 

The necessary condition for solving the problem is to 
derive information about: (i) the probabilities of differ¬ 
ent configurations and (ii) the Maxwell-Boltzmann equi¬ 
librium for classical particles, which is the case that I will 
solve. This approach was already considered by other au¬ 
thors for the largest scales in the Universe (e.g. Saslaw & 
Hamilton 1984; Saslaw 1985) and I also believe that this 
condition could be applicable in some gravitational sys¬ 
tems and in some scales (not all though, see Betancort- 
Rijo 1988): for example, the distribution of galaxies may 
be a good candidate to consider relaxed in which we will 
develop an example of the expression derived here in Sect. 
VIII. Nevertheless, there might be other gravitational sys¬ 
tems for which this approach may be valid. 
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2. The probability of mass and position for a par¬ 
ticle with a known mass function distribution 
and Boltzmann equilibrium. 

Hypothesis 1: We have a grand canonical ensemble 
of N point-particles, interchangeable with an external 
reservoir, with positions ri,...,r jq momenta pi,..., pn 
and masses mi,...,mN respectively. This means that 
the distribution of positions and momenta in the par¬ 
ticles follows a Maxwell-Boltzmann distribution, i.e. 
they are relaxed. 

The probability of a configuration in positions and 
masses is 'PAr(ri,r^r; toi, mAr). First, we define a 
general distribution of masses, and then, once we have 
the mass of each particle established a priori, the prob¬ 
ability of a particle occupying a position depends on its 
mass: 


O 

where TL is the Hamiltonian of the gravitational system, 
U the potential energy and \ ^ kinetic energy. 

We integrate over the momentum space Pi, thus obtaining 
the probability for the positions space configuration: 

X y (3) 

The masses rrii and (3 are constant so the integrals in 
(||) are constant and 


= e^'"^PAr(ri, ...,rN/mi,.. .m]sr)^N{jni,...,m]sr), (1) 

where p is the chemical potential, /9 is a constant, is 
the mass distribution function, and Pn the position distri¬ 
bution function which also depends on the masses as pa¬ 
rameters (the probability of a configuration of positions is 
conditioned by the distribution of masses; the slash stands 
for “conditioned”). We have decomposed the probability 
Vn into the product of three probabilities: first, the prob¬ 
ability of having N particles in a grand canonical ensemble 
(see for example Saslaw 1985, ch. 34); secondly, the prob¬ 
ability of a mass configuration, and, finally, once we 
know the mass of each particle, the probability of the po¬ 
sitions configuration. In the following subsections I obtain 
Pn and $jv. 

2.1. Probability of positions 

We have assumed that the number of particles is 
discrete and that the probability of that they be 
configured in a certain way is given by a function 
PNiri, ...,rjv/TOi, ...,mN). 

To analyse the form of Pn an assumption is needed on 
how position space (or the phase space that includes it) 
is populated, and this can be carried out by assuming a 
state of thermodynamic equilibrium, which was our first 
hypothesis. The Boltzmann distribution is characterized 
by a configuration probability proportional to the negative 
exponential of the Hamiltonian of the system, and the 
Hamiltonian is the sum of two terms, one depending on 
the momenta of the particles and other on the positions 
of the particles. 

The Boltzmann probability distribution function for a 
fixed time (see for example Stanley 1971) is 


2.2. Probability of masses 

Obviously, if we have to take into account the difference 
of masses between particles, we must first derive the mass 
distribution. The usual tool for doing this is the mass dis¬ 
tribution function for a particle, (j}{m), i.e. the probability 
of a particle with mass between m and m-\-dm is (j)(rn)dm. 
Actually, (p represents an average over space of positions 
of the mass distribution. We could consider the distribu¬ 
tion of mass to vary at different positions, i.e. the centre 
of a galaxy cluster with respect to another position, an 
effect which should indeed be taken into account. But this 
is wrong. Remember that by means of (|^), we first chose 
the masses of the particles according to their probabilities, 
and then we assigned a position to the particle conditional 
on having a previously established mass. 

So the probability of N particles having masses 
mi,...,mN, respectively, is 

^N{'mi,...,mN)cc(j){mi)...(p{mN). (5) 


2.3. Normalized probability 

Expressions (|^), (^ and (^), once normalized, lead to 
VNimi, ...,mN;U{ri, ...,rjv/mi, ...,tojv)) 


= e^^^^ (j){mi)...(t>{mN) --, ( 6 ) 

and the partition function of this canonical mass- 
dependent distribution is given by 






dmi... / druN / dr i...dr at 


N=0 


X(j)(rni)...4>{mN)e 


PN{ri, ...,rjv;Pi, ....Pat/toi, ...,mN) 




( 7 ) 
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3. Correlation functions 

The mathematical tools that allow important filtered in¬ 
formation to be extracted from the distribution of parti¬ 
cles are the correlation functions. These will enable us to 
ascertain whether the distribution is Poissonian, or how 
large the difference is with respect to a Poissonian distri¬ 
bution. 

Hereafter, I define certain functions that contain the 
symbol < (• • •) > as meaning an average of the quantity 
(•••). There are two different ways of calculating these 
averages: (i) by volume, and (ii) statistically-mechanically. 
I now describe the two ways: 

(i) In this case, the average is calculated according to 
the expression 


g-/3C/(ri,...rjv/mi,...'mjv) 

x4>{mi)...(l>{mN) - - -. (10) 

In order to simplify the notation, the mass dependence 
of U will not be indicated in what follows. A straightfor¬ 
ward calculation shows that 

< n{r)n(r') >= ^ ^ 

If we assume isotropy, the function only depends on |r — 
r'l; I denote this as < nn > (|r — r'\). 


<(---)>=^ (8) 

where (• ■ •) represents the quantity whose average we wish 
to obtain. It is clear that we are calculating a volume 
average. This is the method that describes the distribu¬ 
tion macroscopically. The observational correlations are 
extracted with this algorithm. 

(ii) In this case, the average calculation is conducted 
by means of the expression 

OO 

<( ■■)>=E 

iV=0 


J dmi...dmN J dri...driq[- 


3.2. Mass-mass two-point correlation function 

The mass density is p{r) = — fi)- The two- 

point correlation function for mass without autocorrela¬ 
tion is 

CXD 

< P{r)p{r') >= 

xp{r)p{r')<l){mi)...4>{mN) -^- 


dmi... / dmAT / dri...dr 


N 


N=0 




dmi... / dmAT / dri...drM 


xVN{ri,...,rN]rni,...,mN), (9) 

where (• • ■) is the same as in (i). This method appeals to 
the microscopic properties of the physical system. 

The definitions that follow include an expansion with 
the second averaging method. 


N 


X 'Y - ri) 5 {r'- rj)^ (l){mi)...(j){mM) 

„-/3C/(ri,...,rjv) 


= g ^0un N{n-i) rf 

N=0 •' 


3.1. Particle-particle two-point correlation function 

The density of the particles is n{r) = ^(’^ ~ ’’i)- 

The two-point correlation function for particles with¬ 
out autocorrelation {i j in the following expression) is, 
from (||) and (^, 

OO 

< nir)nir') >= Y^ e^^^ 

7V=0 


X J dr3...dr7vmim2^i(mi)...^(mAr)e“^*^^’"’’"’’'^’■■■’’"'^^(12) 

With isotropy, the function depends on \r — r'\. 

3.3. Mass-particle two-point correlation function 

With the same philosophy in mind, we define other func¬ 
tions as the two-point correlation function mass-particle: 


dmi... / dmN / dri...dr 


N 


xn{r)n{r')(j){mi)...(j){m]sf) --- 


< p{r)n{r') >= Y 




dmi... / dm^ 


iv=o 


dr 1 ...dr at 


N^O 


^0fiN / / dri...drAr 


N 


Y S{r-ri)S{r' 



N 


X I Y^ miS{r — ri)6{r' 



X(j){mi)...4>{mN) 


g-0U{ri,...,rN) 


0uN^{^ - 1) 


X 
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X / dmi... / dmN / dr3...drN'mi(j){mi)...(j){mN) 


By this hypothesis, the total potential energy is the 
sum of the potential energies between pairs of particles 
{U{ri, ...,rN) = T,^<j - Vj)). Hence, 


xe 


—/3C/(r,r',r3,...,rAr) 


(13) 


Vri < n(ri)u(r2) >= ^ 


„l3nN —PN{N — 1) 


With isotropy, the function depends on \r — r'\. 


N=0 


Z 


dmi... 


3.4- Mass-mass-particle three-point correlation function 
This function is defined as follows: 

< p{r)p{r')n{r") > 

= ^ e^^^ f dmi... f dmN j dri...dr 


X J dmN J dr3...drjv^(TOi)...(/>(mjv)e dU{ri,...,rN) 


N 


N 


N=0 

( 


N 


mimjS{r — ri)5{r' — rj)6{r" — Vj) 


:,j,k=i 


X(j)(mi)...4>{mN)- 


,-pU{rx,...,rN) 


Vr-1 Vi2(ri - ■^2) + X! - Vi) 

i^3 

This becomes 

V,, < n(ri)u(r2) >= -^iN - 1) j 

X J dmN'^riVi2iri - r2) J dr3...drN(l){mi)...4i{mN) 


(16) 


Ar=0 


= —— J dmi... J dmN j dr^... x j dmi... j dmN J dr 3 V r^Vi^iri - r^) 


xdrNmim24>imi)...(l){mN)e 


—/3U{r,r',r",r4,...,rN) 


(14) 


With isotropy, this function depends on three variables: 
\r — r'\, |r — r"| and |r' — r"| . 


4. A force equation for homogeneous, isotropic 
distributions in Boltzmann equilibrium 

In this section, I derive the force equation and then intro¬ 
duce a Newtonian form for the interaction. 

The only trick is to apply the operator over ex¬ 
pression & and to develop the expression 

Vr, < n{ri)n{r2) >= ^ f 

N=o '' 

X J dmN J dr3...drN4>{mi)...(j){mN)e~^^^''^'"'’'^^^ 


X J drii...drN4’(mi)...(j)(mN)e 

Hypothesis 3: We assume homogeneity and isotropy. 

This means that < n{ri)n{rj) > only depends on distance 
rij = \ri — rj \ and V{r) = V{r). The operator can be 
expanded and thus (I note < nin 2 > instead of < nu > 
to avoid confusion; the same with < pp > and < ppn > 
in the following equations) 


d < nin2 > (ri2) 

M2-6' 

N=0 


dri2 


z 


X J dmi... J dmN[ri 2 ^^^^^^-^] J dr 3 ...dr n 


X(j){mi)...4>{mN)e dH(ri,...,rjv) _|_ ^ 
-I3N{N -1){N -2) 


N=0 


dr 3 / dmi... / dmN 


x{-(3)VriU{ri, ...,rN)- 


(15) 


Hypothesis 2: The interaction among particles is be¬ 
tween pairs of particles, and is proportional to their 
masses, as a function of their distance (it is a central 
force). 


ri2 cos(ri2,ri3) 


dVisjris) 

dri3 


X J dr4...dr7v(/)(mi)...(/)(m7v)e dH(ri,...,rw)^ 


(18) 
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I rename r = ri 2 , s = ri 3 , 9 = angle between ri 2 and 
ri 3 and assume that Vij{r) = mimjv{r) (valid in the light 
of Hypothesis 2) and substitute some expressions using 
and taking into account the isotropy assumed 
in Hypothesis 3: 


— 1 5 < 711712 > (r) 
[3 dr 


dv{r) 

dr 


< PP> (r) 


+27r 



dO sin 6 cos 9 



ds 


dv(s^ / - 

X— < pip 2 n 3 > (r, s, V — 2rs cos0). (19) 

This is the force equation. It gives us a relationship be¬ 
tween different correlation functions in a distribution that 
follows our hypothesis and the force that is represented by 
means of v. 


5. Newtonian interaction and expansion 
5.1. Newtonian interaction 

Before proceeding, I will comment on some aspects regard¬ 
ing the Newtonian case in particular. Certain problems are 
associated with the application of principles such as ther¬ 
modynamics and statistical mechanics in groups of parti¬ 
cles with r“^-type force (Taff 1985): some divergences are 
found, and there is non-saturation of gravitational forces 
(Levy-Leblond 1969). Some authors take the view that 
there can be no rigorous basis for applying statistical me¬ 
chanics in such a system (Fisher & Ruelle 1966). This re¬ 
sult is not a consequence of the r“^-type force but rather 
of its unshielded character (see Dyson & Lenard 1967 for 
a discussion of the electrostatic case). 

In any event, one avoids the non-self-consistency of the 
problem by truncating the integral limits at a hnite radius, 
or by cancelling the correlation functions that fall under a 
given lower unit, assuming points with negligible volume 
as is the case in real physical problems. In my opinion, 
A^-body systems exist in nature under Newtonian gravi¬ 
tational forces, and to invoke a distribution of particles in 
such systems is not necessarily inconsistent. The possible 
infinities that appear in some expressions are only math¬ 
ematical problems which are not present in nature and 
can be solved once we give our data the conditions that 
will accommodate our physical reality to a mathematical 
model (I mean to avoid the infinite proximity of particles 
by means of a cut-off, etc.). We can look at it from an¬ 
other point of view: a statistical thermodynamical system 
in equilibrium cannot achieve singular states with infini¬ 
ties, except as a set of zero measurements, because the 
probability of getting a singular state, with a very small 
distance between some particles, would take an extremely 
long time, as if approaching a singularity. 


Expression (H) has been obtained regardless of the 
force type, so it possesses a general validity. Now, when 
we introduce a Newtonian gravitational force the expres¬ 
sion continues to be valid. If an infinity appears in the 
next expression it is only a question of truncating the in¬ 
tegrals or selecting the best correlation function that does 
not produce divergences. When we take the cut-off, we 
neglect the probabilities near the singularity, a set with 
dimensions greater than zero but small enough. The intro¬ 
duction of cut-offs will make the results dependent upon 
the details of the regularization, so the selection of the 
cut-offs must have a physical basis. 

Hypothesis 4: The interaction between pairs of particles 
is the Newtonian gravity force. 

This hypothesis obviously includes Hypothesis 2. With a 
Newtonian potential: Vij{rk — ri) = Then, 


dv{r) _ G 
dr 


We set C = iPG) ^. This leads to 


d < 771772 > (r) 1 

^ - = ^12 < > (^) 


-|-27r / sin 6* cos 6* / ds 

Jo Jo 


X < PlP2I^3 > ^/r"^ + s"^ — 2rs cos 9) , (20) 

which gives us a relationship among certain correlations 
of the distribution and C. 

The proof that the system is valid for achieving ther¬ 
modynamic equilibrium can be found in Lieb & Lebowitz 
(1973), where a general Coulombian system is considered. 
Nevertheless, this does not imply that all Newtonian grav¬ 
itational systems are in equilibrium. 

5.2. Effects of the expansion of the Universe 

When we consider the galaxy distribution in the large- 
scale structure of the Universe, we must bear in mind the 
expansion, so: 

Hypothesis 5: The system of particles is distributed over 
a space in expansion. 

We must use comoving coordinates in order to maintain 
the zero peculiar velocity as the most probable one (be¬ 
cause the above formulation gives the probability as pro¬ 
portional to e“^, where p is the momentum). We consider 
the proper motions, not the background: the kinetic en¬ 
ergy derived from peculiar velocities and the comoving 
potential energy. 

It is derived in Saslaw & Fang (1996) that considera¬ 
tion of the expansion is equivalent to taking into account 
the gravitational effects of the local fluctuating part of the 
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density field when we use comoving coordinates, i.e. we 
should substract the mean density from the density held. 
Since the mean density does not depend on r or 0 and the 
integral dd sin 9 cos 9 = 0, equation ( |20| ) is modihed by 
simply replacing < pip 2 > (r) by < pip 2 > (r)— < p >^: 


d < nin2 > (r) 1 2^ 

-C - - -^ = ^(< P 1 P 2 > (r)- < p >^) 


-\-2tt f d9sin9cos9 
Jo 



ds 


X < PiP 2 ns > (r, s, \/r2 -I- - 2rs cos 9). 


( 21 ) 


6. The meaning of (3 

We know the meaning of (3 in statistical physics; it is equal 
to where Kb is Boltzmann’s constant, and T is the 
absolute temperature. In the large-scale (or any other) 
structure we can establish a similar meaning by developing 
a kinetic theory, in which the particles are the galaxies 
themselves (or other particles). 

If we integrate expression (H) over position space, we 
get 


7. How to obtain the correlations from the distri¬ 
bution 

To this end, we have to use eq. (||) as dehned previously. 
However, when we have a discrete number of points in¬ 
stead of a continuum distribution, perhaps it might be 
better to use other equivalent expressions. 

When homogeneity is given, one method discussed by 
Rivolo (1986) is to use the estimator 


< niU 2 > (r) = 


<. n > Ni (r) 

N ^ VAr) ’ 

i—1 ^ ' 


(26) 


where Ni{r) is the number of particles lying in a shell of 
thickness Sr from the ith particle, Vi(r) is the volume of 
the shell lying within the sample volume and < n > is the 
average density in a macroscopically homogeneous system 
{< n > is independent of position). 

The evaluation of < pip 2 > and < pn > would be: 


< PiP 2 > (r) 


< p> miMi{r) 

M ^ U,(r) 

1—1 ^ ' 


and 


(27) 


V{pi, ...,pAr;mi, ...tojv) 

(X . ( 22 ) 

The mean value of the peculiar velocities for all clus¬ 
tering galaxies is 

_ f dm(j){m) J dp^e~^^^ 
f dm4>{m) f dpe~i^^ 


J _ / 8 \ f dm(j){m)m 

We can also express this as 


V = 


irp J 


(24) 


where p = ( ^ (j represent the same 

constant as in the force equation for the Newtonian case. 
Hence, note that the parameter G, which we could obtain 
in the equality will give information about the ve¬ 
locity held of the particles, or, vice-versa, we could obtain 
G for (|^) from the mean velocity. 

The parameter (3 is even more directly related to the 
kinetic energy. An analogous derivation leads to 

3N 8 

Euin = ^ = (25) 


< Pi «2 > (r) 


<n> miNi{r) 

M ^ mr) 
1—1 ^ ' 


(28) 


where M = Mi{r) is the mass lying in a 

shell of thickness Sr from the zth particle. Knowledge of 
rrii is sometimes problematic and its solution differs in 
each case. In the case of galaxies or stars information is 
required about the mass-luminosity relationship as well as 
data on magnitudes and distances. This could also be pro¬ 
vided with a knowledge of the mass distribution function 
according to different zones, ( (j){m) would be the average 
of the different mass distribution functions in the entire 
space) and by randomly assigning a mass rrii to each par¬ 
ticle following these distribution functions that would give 
us not the real mass distribution but an equivalent one. 
In any case, the problem of assigning mass is different in 
each case and generally requires tailor-made solutions for 
each one. 

The three-point correlation function is also obtainable 
by counting groups of three particles with different dis¬ 
tances between them. Usually, in the isotropic case, this is 
approximated as a function of different two-point correla¬ 
tion functions. For example, in liquid theory the so-called 
superposition approximation is commonplace (see for ex¬ 
ample March & Tosi 1976 and first formulation of it in 
Kirkwood 1935), and applied to the mass-mass-particle 
case would take the form 


< PiP2n3 > {r,s,t) 


< P 1 P 2 > (r) < Pi712 > (s) < P2n3 > jt) 
< p >2< n > 


a fundamental result of statistical mechanics. 


(29) 
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From the large-scale distribution of galaxies in the Uni¬ 
verse, we can extract statistical information (as in Saun¬ 
ders et al. 1991). When we observe the projected distri¬ 
bution onto a 2-dimensional surface, i.e. we do not know 
the distance of the objects, we can obtain the correlation 
on the 2-dimensional surface (angular correlation) and re¬ 
late it to the 3-dimensional distribution correlations by 
means of Limber’s equation (Peebles 1980). Maddox et al. 
(1990) obtained the two-point angular correlation func¬ 
tion on large scales, so we can derive < nin 2 >■ To assign 
the masses in order to achieve < pip 2 > and < pin 2 > is 
again troublesome and depends on the data available on 
the masses of the galaxies. 

For the mass of galaxies in the large-scale structure, fu¬ 
ture data will become available with the help of the DENIS 
project (Mamon 1995), a near-infrared sky survey. The 
near-infrared is thought to be a better tracer of the stel¬ 
lar mass in the galaxies, and if the stellar mass content 
follows the total mass content (including dark matter), 
near-infrared surveys should be the best way of obtaining 
the distribution of matter in the Universe. 

Also, on other astronomical scales we find point- 
particle distributions and we may obtain their correla¬ 
tions. Borgani et al. (1991) obtained the correlation func¬ 
tions for scales between 3 and 350 kpc, and this could be 
done for smaller scales as well. Another problem is posed 
by considering our hypotheses as valid, especially Hypoth¬ 
esis 1. 

An application to any other distribution is also possi¬ 
ble. 


8. Application on the Large-scale distribution of 

galaxies 

The actual application of these tools is quite difficult since 
real data are not easily available. To illustrate the way to 
proceed with this method, we apply the above theoretical 
results to a practical case with observational data, in the 
A-body system of the large scale distribution of galaxies 
in the Universe (see for example Peebles 1980; Borgani 
1995) which is a homogeneous and isotropic distribution. 
A recent model, using equilibrium statistical mechanics as 
well as other considerations, was also developed in Perez- 
Mercader et al. (1996). 

We will assume the validity of the assumptions made 
in this paper and adopt an additional assumption here 
for the particular case of the use of (^), the superposi¬ 
tion approximation, and biasing of the p fluctuation as 
proportional to the n fluctuation where the constant of 
proportionality depends on the scale: 


6n 6p 

— = b{r) — , 
n p 


(30) 


where b is the biasing parameter, dependent a priori on r, 
the distance between both particles. Astronomers usually 
make a stronger assumption by taking b to be constant. 


but I am not going to be so restrictive. I adopt a unit 
system where < p >=< n > and set < nin 2 > (r) =< 
n (1 + C{r)). As a consequence, 

< Pin 2 >=< n >'^ {1 + b~^^{r)) (31) 

and 


< PiP2 >=< n (1 -b 5 
Thus, (iB leads to: 


-Ci'ir) 


j,2 


-\-27t < n > 



ds[l-|-6 ^(s)^(s)-|-5 ‘^{r)^{r) 


+b ^(r)6 ^ {s)^{r)^{s)] f dO sinO cos9 

Jo 


(32) 


xb ^(\/— 2rscos0)^(\/-|- _ 2rscos6). (33) 

This expression resembles the one given by Peebles (1976) 
for 6=1, the cosmic virial theorem, but with further 
generality because expression (^) takes into account that 
< nin 2 > (r) piP 2 > (r). 

An observational ^ was achieved by Groth & Peebles 
(1977): 


C(0 


^^ for 0.36 ^ Mpc < r < 106 ^ Mpc 


(34) 


where ro = 4.16“^ Mpc and 7 = 1.77. The parameter ^ is 
negligible for r > 106“^ Mpc, so it is taken as zero. Also 
a cut-off is taken for r < 0.36“^ Mpc because separated 
galaxies cannot be at distances less than this cut-off (if 
two galaxies have a distance between themselves less than 
this then they are considered as only one galaxy), and this 
ensures convergence of the integral in ( |^ ) . 

In order to estimate the behaviour of 6(r), we introduce 
it as a power-law dependence, such that 


b{r) = j (35) 

and ht the best values of k and r* to solve (|^. The 
constants are derived from observational data: < n >= 
0.02 Mpc“^ (Allen 1973), 6 = 0.60 (the actual value of 
the Hubble constant in units of 100 km s“^ Mpc“^) and 
V ~ 100 km s“^ (Allen 1973). 

If we now wish to obtain C, according to (p^, we must 
convert G from our units (time unit: second; length unit: 
1 6“^ Mpc; mass unit: the average mass of a galaxy, Wgai 
in K‘kg) to the MKS system, multiplying by a factor of 
(3.08 X 10^^)“^TOgai. Thus, denoting by Gmks the value 
of G in the MKS system, we have 


V = 


^ 8Gmks(3.08 X 10"^)-3mgaiG 


1/2 


6 ^ Mpc s ^,(36) 


TT 
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so with Gmks = 6.672 x 10^^, Mq = 
1 Mpc = 3.08 X IQi® km 

1.99 X 103O kg and 

_ 3.45 X 1013 Mq 
o ~ 

m 

(37) 

This, together with an estimate of m » 
1973), gives us 

^ 8 X 10i° Mq (Allen 

C Ri 4.3 X 10^. 

(38) 


The introduction of all these data into (^) and some 
calculations lead to the following expression for the al¬ 
lowed r values: 


n on „ in3 -2.77 , n i 2^ -3.77-2 k 0.291r-J 

9.25 X 10 r — 12 .r — - 


r+10 


ds, „ 

— [1 + 12.1rJ 


-1 

-l.77-n 


s < 0.3 
0.3 < s < 10 
s > 10 


k= 0.52, r =8.0 



+ 12.Ir^ r 


' +83.1r:;: 


-1 

-1.77-k 


s < 0.3 
0.3 < s < 10 
s > 10 


Fig. 2. Comparison between the left- and right-hand side of 
expression (p^ for best fit of k and r*. 


2.23 - 


7 ({ 


0.3 

|r- 


2.23-« 
2.23-k 


|r — s| < 0.3 
|r — s| > 0.3 


- (r + s) 


2.23-7 


0.23 • 


({ 


0 3u.^o- 

|r - s|° 


|r - s| < 0.3 \ 
jr - s| > 0.3 / 


- (r -I- s 


,0.23-^ 


K=0 (without biasing) 


The best fit, calculated numerically, is shown in Fig. 
2 , where left- and right-hand sides of expression for 
parameters k = 0.52 and r* = 8.0 h~^ Mpc are plotted. 
( 39 ) With these parameters we get 

/ \ 0.52 

• <«> 

or, using a different expression 



r (h ’Mpc) 


Fig. 1. Comparison between the left- and right-hand side of 
expression (M) for the non-biasing case. 

If non-biasing, i.e. 6=1, were assumed for all scales, 
we would obtain the results plotted in Fig. 1, so biasing is 
neccesary. 


/ \ “ ° 

= ( safe-l MpJ («) 

in the range above defined for where < pp > (r) =< 
n {1 + ^pp{r)). The obtained power law for mass cor¬ 
relation, “the —2.8 power law” is different from that of 
objects’ correlation, “the —1.8 power law”. Of course, the 
left and hand sides of were not expected to agree 
perfectly and in fact they do not (Fig. 2) because we as¬ 
sumed a 6 dependence, a power law, which might not be 
very realistic. A worthwhile result is that a deviation from 
6 = 1 is necessary to get an agreement between the left- 
and right-hand sides of (^), and the solution must be 
close to 6(r) = 0.34r°-^^. This result increases the informa¬ 
tion about the mass-mass correlation function (see Peebles 
1980), which was unknown till now. If this result is true 
and is verified by alternative methods, it would give us 
important information about the “dynamics” in the large 
scale distribution of matter in the Universe. 

The dependence of the outcome on the cut-off value 
is not negligible, but neither is it very pronounced. Some 
numerical results were obtained with other values of the 
cut-off, and the qualitative result does not differ too much: 
a cut-off at r = 0.1 h~^ Mpc instead of r = 0.3 h~^ Mpc, 
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gives us K = 0.47 instead of 0.52, and r* = 8.1 h ^ Mpc 
instead of r* = 8.0 h~^ Mpc. 



r (h ' Mpc) 


Fig. 3. The galaxy-galaxy correlation function, derived from 
observational data and ^pp, mass-mass correlation function, de¬ 
rived from ^ through the force equation explained in this paper. 

The direct consequence of this is that mass is more 
correlated than the objects (Fig. 3). This means that dark 
matter must exist near galaxies and clusters of galaxies to 
increase the density contrast for short distances from an 
object (dark matter in form of discrete unseen galaxies, for 
instance dwarfs galaxies), unless the most massive galaxies 
are clumped together or there is some other solution, but 
something must explain the difference in both correlation 
functions for galaxies and mass. Deviations from equilib¬ 
rium might also be responsible for part of the deviation 
attributed to biasing, as far as it is only an approximation. 
I am unware of how this may affect my results, essentially 
because the extent of this deviation is unknown, though 
I expect it to be not very large. A large departure from 
equilibrium would produce a fast evolution of the distri¬ 
bution and the distributions of matter far from us would 
have different distributians. This latter is not observed for 
long-range correlations u, so we must infer that the depar¬ 
ture from equilibrium conditions cannot be too high unless 
perhaps for very small ranges. For typical scales in galaxy 
clusters, we also find strong evidence of equilibrium (Carl- 
berg et al. 1997). This discussion is beyond the scope of 
this paper. In Saslaw & Hamilton (1984, their Sect. 6), we 
find further arguments in favour of this: “Gravitational 

^ Broadhurst et al. (1990) deduces from observational data a 
regular distribution of galaxies similar to the nearby structnre 
up to distances of ~ 1000 h~^ Mpc (a back in time further 
away from ^ 4 x 10® yr) 


iV-body experiments give good agreement with the the¬ 
ory. This shows that even though an equilibrium theory 
may not explain the ultimate fate of galaxy clustering, it 
does provide a good description over extremely long time 
scales when the correlations are ‘frozen out’ by the expan¬ 
sion of the Uni verse... [T] he equilibrium theory seems to 
explain most of the results.” 

The main question is how equilibrium could be reached 
in a short lifetime of the Universe. Violent relaxation 
which enormously decreases the relaxation time is a pos¬ 
sible solution. Indeed, Saslaw (1985) points out that this 
must be the mechanism that governs the system due 
to large-scale collective modes (see chapter 38 of Saslaw 
1985) and Henriksen & Widrow (1997) makes numerical 
simulations achieving this. The consideration of a steady 
state in the large-scale structure or a relaxation included 
in the initial conditions of the large-scale structure dynam¬ 
ics, before the formation of the galaxies, are other possible 
explanations. 

The deviation from equilibrium and the corrections to 
make to our equations to take these effects into account 
are topics for future papers. Further research is necessary 
in this area to render these results more accurate. With 
this example, we wanted to show the way of working with 
the expressions described in this paper. 

Readers might ask why I have not used the force equa¬ 
tion to derive the two point-correlation function. To do 
this, I would need a knowledge of the biasing first, I can¬ 
not derive both things at the same time. Since the two 
point correlation function is better known that the bias, 
I decided to apply the method as described above. I had 
previously done some calculations to derive the two-point 
correlation function assuming non-biasing and the result 
was not compatible with observations, so I rejected this 
hypothesis thinking that biasing is necessary as I obtained 
here. 

9. Results, other applications and further com¬ 
mentaries 

We have a numerical relationship between the distribution 
of galaxies in space, which is represented by the correla¬ 
tion functions, and the two-galaxy interaction, which is 
represented by the potential energy. 

We have obtained the equality (iil) that must be fol¬ 
lowed by the distributions under Hypotheses 1 to 5. Also, 
we have equality (H) for the correlation functions and any 
interaction force under Hypotheses 1, 2 and 3 (we could 
also obtain an expresion like this including the expansion 
of space by means of the method explained in the subsec¬ 
tion dedicated to the expansion). 

Equation (|^ ) relates the distribution correlation with 
the mean velocity of the galaxies by means of C with (Q), 
and I believe it will be useful for obtaining a parameter 
from others that are already known in the distribution: the 
mean velocity from a complete knowledge of the correla- 
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tion functions, or an unknown parameter in the correlation 
function from the rest of the data. We could even obtain 
more than one parameter: two or three (in my opinion, 
more than three are too many) that follow the equality 
between the right- and left-hand sides of equation ( pl|). 

When equality between the two-sides of equation (^l[) 
is unattainable, this will indicate that our hypotheses are 
unsuitable. Probably, the most doubtful hypothesis is the 
first, i.e. that of Boltzmann equilibrium, ancLit is possible 
to verify the relaxation using this equation. El 

In a sufficiently evolved system, Boltzmann equilib¬ 
rium is achieved because the particle-points are classical 
particles and the probability of a state in such a case is 
proportional to the number of different states for each par¬ 
ticle that preserves the number of particles and the total 
energy (the reader is referred to any book dealing with the 
foundations of statistical mechanics, e.g. Tolman 1938). It 
is also true that after a long time, the systems become 
virialized, and many systems are known to be in these 
conditions, although not all of them. 

Otherwise, equation the force equation^ is a tool 
for looking for the kind of interaction in a system follow¬ 
ing Hypotheses 1, 2 and 3 (as purposed in Goldman et al. 
1992). Once we have derived all the parameters of the dis¬ 
tribution we can fit a shape for U, a two-body interaction 
potential that fits the equality, with or without expansion. 
We could even obtain other unknown parameters (C for 
example). 

In order to demonstrate what might be the caveats in 
the implementation of the force equation, I developed a 
real example in the previous section where it was used to 
infer information about the mass-mass correlation func¬ 
tion from the galaxy-galaxy correlation function, and the 
average density and peculiar velocity in the large-scale 
distribution of galaxies in the Universe. Further improve¬ 
ments are necessary, both in the observations and the the¬ 
oretical assumptions, to obtain an accurate result, but the 
method is at least capable of telling us that the mass is 
more correlated than the galaxies at short distances (Fig. 
3) when we assume relaxation on scales greater than 1 h~^ 
Mpc. 
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